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IGNITION,  COMBPSTION,  DETONATION,  AND  QUENCHING 
OF  REACTIVE  GAS  MIXTURES 


(Five  Different  Topics  of  this  Research  Program  Have  Been  Investi- 
gated During  this  Report  Period  — 1 April  1975  to  3I  March  1976) 


I.  EFFECT  OF  INITIAL  TEMPERAIURE  ON  THE  DETONATION  INDUCTION 
DISTANCES  OF  HYDROGEN-AIR  MIXIURES 


Although  many  papers  have  been  published  on  the  structure  of  caie- 
dimensional  detonation  waves  propagating  into  combustible  gas  mixtures 
at  room  temperature,  hardly  any  studies  have  been  made  on  the  effect  of 
initial  temperature  on  the  initiation  process.  Several  years  ago  Edse 
and  Lawrence  Jr.^  measured  the  detonation  induction  distances  in  low- 
temperature  hydrogen-oxygen  mixtures.  The  present  study  was  undertaken 
to  determine  whether  fuel-air  mixtures  behave  like  fuel-oxygen  mixtures. 


A.  APPARATUS  AND  MEASUREMENTS 

A 6.4  m long  stainless  steel  tube  with  an  inner  diameter  of  5 cm 
and  a weill  thickness  of  0.635  cm  was  used  for  these  exx>eriments . Fig- 
vire  1 is  a photograph  of  the  tube,  and  the  locajsions  of  the  probes  are 
given  in  Fig.  2.  The  detonation  tube  was  encased  in  a 0.635  cm  steel 
^Jacket  which  was  welded  to  the  tube  and  insulated  to  the  surroundings  by 
two  layers  of  Armstrong -Armaflex  rubber,  each  having  a thickness  of  I.9  cm. 
A Minnesota  Valley  Engineering  (mvE)  Model  VGL-I6OL  liquid  nitrogen  con- 
tainer was  used  to  pass  liquid  nitrogen  through  the  jacket  at  a speed 
which  established  the  desired  low  temperature  of  the  combustible  gas  mix- 
ture in  the  detonation  tube.  This  mixture  contained  30^^  hydrogen  and  70^ 
dry  air  by  volume.  Complete  mixing  was  established  in  a specially  de- 
signed mixing  chamber  from  which  the  mixture  was  transferred  directly  as 
a continuous  flow  to  the  detonation  tube.  Prior  to  this  transfer  the 
detonation  tube  was  evacuated  to  remove  all  moistvire. 

The  initial  temperatvire  of  the  combustibe  gas  mixture  was  measured 
by  two  copper-constantan  thermocouples.  A 0.I3  mm  pyrofuse  wire,  ex- 
ploded by  a 40  V ac  power  supply,  was  used  to  ignite  the  hydrogen-air 
mlxt\ire.  Two  Kistler  pressure  transducers  were  employed,  together  with 
a Tektronix  lype  555  dual  beam  oscilloscope,  to  determine  the  detonation 
wave  speed  and  to  measure  the  pressure  at  two  locations  in  each  experi- 
ment. The  transducers  were  placed  at  various  distances  from  the  ignitor. 

At  least  three  measurements  were  made  for  each  location.  Three  typical 
oscilloscope  traces  are  shown  in  Figs.  3®~c.  All  experimental  data  are 
tabulated  in  Tables  I,  II,  and  III  and  depicted  graphically  in  Figs.  4 
and  5»  The  detonation  induction  distances  as  obtained  from  Fig.  5 are 
compiled  in  Table  IV. 
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3(c)  Initial  Tempera- 
ture = 223  K 

Upper  Trace;  I5  cm 
from  ignitor  sweep 
rate  = 5 ms/div  ver- 
tical sensitivity  = 
4.6  atm/div 

Lower  Trace;  45  cm 
from  ignitor  sweep 
rate  = 5 ms/div  ver- 
tical sensitivity  = 
5.0  atm/div 


Fig.  3 - Typical  pressure  traces  in  the  detonation  induction  region  at 
various  initial  temperatures 


3(a)  Initial  Tempera- 
ture = 123  K 


Upper  Trace;  15  cm 
from  ignitor  sweep 
rate  = 5 ms/div  ver- 
tical sensitivity  = 
4.4  atm/div 


Lower  Trace:  45  cm 

from  ignitor  sweep 
rate  = 5 ms/div  ver- 
tical sensitivity  = 
4.2  atm/div 


Lower  Trace;  45  cm 
from  ignitor  sweep 
rate  = 5 ms/div  ver- 
tical sensitivity  = 
6.3  atm/div 


Upper  Trace : 15  cm 

from  ignitor  sweep 
rate  = 5 ms/div  ver- 
tical sensitivity  = 
4.6  atm/div 


3(b)  Initial  Tempera- 
ture = 173  K 


Table  1.  Wave  speeds  and  wave  pressures  in  the  detonation  induction 
region  for  a 30^  hydrogen  in  air  mixture  at  an  initial  tem- 
perature of  123  K and  an  initial  pressure  of  1 atm 


Distance  from 
Ignitor 
(cm) 

Wave  Speed 
(m/s) 

Wave  Pressure 
(atm) 

Remark 

15 

45 

— 

7 ^+2.3 

-2.7 

The  detonation 
wave  observed  at 
242  cm  from  the 
ignitor  decayed 
further  down  the 
tube 

90 

272116 

“•3!o;8 

149 

180 

668±37 

242 

1922115 

18.210.5 

302 

^ ■’-221 

no  y+3*3 

13.7.4  2 

332 

8+282 

^^^-168 

13- 

-2.5 

394 

+222 

12  2*^-^ 
^•^-1.4 

Table  2.  Wave  speeds  and  wave  pressures  in  the  detona- 
tion induction  region  for  a 30^  hydrogen  in 
air  mixture  at  an  initial  temperature  of  173  K 
and  an  initial  pressure  of  1 atm 


Distance  from 
Ignitor 
(cm) 

Wave  Speed 
(m/s) 

Wave  Pressure 
(atm) 

- 

4.610.2 

45 

- 

6.2+° -3 
-0.5 

90 

200+22 

6.510.4 

149 

564:^3 

8.010.2 

180 

242 

302 

332 

1992!®7 

394 

2000+0 

-ft  _+0.4 
l®'°-0.3 

454 

1983+0 

ift  n+°-^ 

577 

197210 

-ft  -+0.4 
l®-°-0.3 

630 

198110 

i8.0_o.3 

- Not  measured 


Graphic  presentation  of  wave  speed  data  from  Tables  2,  and  3 


Table  4.  Average  induction  distances  (dind)  ® hydrogen  in  air 
mixture  in  a 6.4  m long,  5 cm  I.D.  tube  at  various  initial 
temperatures  (Ti):  Initial  Pressure,  Pi  = 1 atm 
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123 

173 

223 

295 

230 

332 

450 

No  Detona- 
tion was 
observed 

Note:  Induction  distance  (dind)  has  been  defined  as  the  distance  be- 

tween the  source  of  ignition  and  that  point  in  the  tube  where 
the  wave  first  travels  at  the  speed  of  the  C.J.  detonation  wave. 
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1.  Determination  of  Normal  Burning  Speeds 

In  order  to  assess  the  influence  of  the  normal  flame  propagation 
rate  on  the  transition  from  deflagration  to  detonation  the  flame  propa- 
gation rates  of  hydrogen-air  mixtures  were  measured  at  various  tempera- 
tures. By  cooling  and  heating  i^he  gas  mixture  it  was  possible  to  make 
measurements  over  a temperature  range  from  123  to  784  K. 

At  temperatures  above  room  temperature  the  straight  burner  tube  as 
well  as  the  individual  gas  flows  were  heated.  The  hot  gases  were  intro- 
duced into  a mixing  chamber  which  was  connected  to  the  heated  70  cm  long 
stainless  steel  burner  tube  having  an  internal  diameter  of  0.458  cm.  To 
avoid  flash-back  a hydrogen-rich  mixture  (53*6?&  hydrogen  and  46.4%  air) 
had  to  be  used.  Photographs  of  the  flames  were  taken  at  original  size 
on  Kodak  Plus-X  35  film.  The  temperature  of  the  unburned  (heated) 
hydrogen-air  mixture  was  measured  at  the  exit  of  the  burner  tube  with  a 
chromel-alumel  thermocouple  just  prior  to  ignition  and  subsequent  expo- 
sure of  the  film. 

The  flame  speeds,  Uf,  given  in  Table  V and  shown  in  Fig.  6,  are 
based  on  the  average  angle  of  the  inner  flame  cone. 

For  the  measurements  at  low  temperatures  the  burner  consisted  of  a 
70  cm  long  copper  tube  with  a 0.635  cm  internal  diameter.  Mixing  was 
acconplished  in  a coiled  copper  tube  with  an  outside  diameter  of  0.635  cm. 
The  diameter  of  the  coils  was  10  cm  and  the  spacing  was  1.3  cm.  In  order 
to  avoid  flash-back  a small  nozzle  with  an  inside  diameter  of  0.3175  cm 
was  silver- soldered  to  the  tip  of  the  burner  tube.  With  this  nozzle  it 
was  possible  to  make  flame  speed  measurements  with  30%  hydrogen  and  70% 
air  mixtvires.  Both  the  bxxrner  tube  and  the  mixing  coil  were  immersed  in 
a large  Dewar  vessel  partially  filled  with  liquid  nitrogen.  The  results 
of  the  measurements  are  given  in  Table  V and  shown  in  Fig.  6. 


B.  DISCUSSION  OF  RESULTS 

Figure  5 shows  the  wave  pressure  in  the  detonation  induction  region 
as  a function  of  distance  from  the  ignitor.  On  the  basis  of  previous 
measurements^  it  has  been  assumed  that  the  maximum  pressvire  attained  in 
the  induction  region  is  practically  the  same  for  all  initial  temperattires . 
Ihese  pressvires  were  also  observed  in  a 42  m long  and  7*5  cm  internal 
diameter  detonation  tube  at  room  temperature  (Fig.  7).  However,  the  tube 
used  for  the  measurements  at  low  temperatvires  was  not  long  enough  to  per- 
mit the  formation  of  a stable  detonation  wave  when  the  mixture  was  ini- 
tially at  room  temperature.  By  ccmparing  the  results  of  Fig.  7 with  those 
shown  in  Fig.  5>  it  can  be  concluded  that  the  pressxire-distance  profile 
within  the  induction  region  at  room  temperattire  is  similar  to  those  at 
low  temperatvires.  The  pressure  overshoot  as  indicated  in  Ref.  1,  which 
is  attributed  to  either  shock  merging^  or  to  the  boundary  layer  type  ini- 
tiaticxi,®  is  also  observed,  before  the  onset  of  detonation.  The  detona- 
tion induction  distances  are  significantly  reduced  as  the  temperature  of 
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Table  5«  Normal  flame  speed  of  hydrogen-air  mixtures  at  various 
initial  temperatures  by  the  flame  cone  method:  Initied. 

Pressure  = 1 atm 


Mixture  Composition 
(by  volume) 

Initial  Temperature 
(K) 

Normal  Flame  Speed 
(m/s) 

123 

1.40 

290 

2.64 

Hydrogen  = 53-6^ 

523 

5.09 

Air  = 46.4^ 

620 

6.24 

755 

11.71 

784 

12.18 

Hydrogen  = 

123 

1.27 

Air  = ICfio 

290 

2.40 

12 


Fig.  7 - Pressure-distance  profile  at  room  temperature 
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the  initial  gas  mixtiire  is  reduced.  At  123  K the  detonation  wave  is 
already  formed  230  cm  away  from  the  ignitor  whereas  no  detonation  was 
observed  at  room  temperature  in  this  tube  (Table  IV). 

Although  at  123  K the  detonation  wave  was  formed  after  230  cm  of 
travel  of  the  deflagration  wave,  it  was  not  stable.  Apparently  heat  los- 
ses to  the  cold  walls  of  the  tube  are  responsible  for  this  decay.  At 
173  and  223  K this  decay  of  the  wave  was  not  observed. 

The  normal  flame  propagation  rates  of  hydrogen-air  mixtures  at 
various  temperatoires  are  compiled  in  Table  V and  depicted  graphically  in 
Fig.  6.  These  data  reveal  very  clearly  that  the  initial  temperature  of 
the  gas  has  a substantial  influence  on  the  b\irning  speed  of  this  gas  mix- 
ture. This  observation  has  also  been  reported  by  Duggar  and  Heimel'*^  for 
various  hydrocarbon-air  mixtures.  The  res\ilts,  as  shown  in  Fig.  6 and 
given  in  Table  V for  hydrogen  and  air  mixtures,  follow  an  empirical  re- 
lationship which  is  quite  similar  to  that  developed  by  Heimel  and  Weast^; 
i.e. , 

Uf.  = B + C*Ti“  , 

where 

Uf  = laminar  flame  speed  (m/s)  , 

Ti  = initial  gas  temperature  (k)  , 

B,  C,  and  n are  empirical  constants,  and 
the  value  of  n is  between  2 and  3« 

In  the  present  case  the  values  of  B and  C are  1.128  and  0.000018,  re- 
spectively, and  the  value  of  n is  taken  to  be  2,  so  that  the  relationship 
between  flame  propagation  rate  and  temperature  for  a 30^  hydrogen  in  air 
mixtvire  becomes 

Uf  = 1.128  + 0.0000l8-Tf  . 

Thus  the  laminar  flame  speed  of  this  mixture  is  approximately  proportional 
to  the  square  of  the  absolute  temperatxire  of  the  initial  gas. 

Ibe  question  that  now  arises  is  "How  can  the  detonation  induction 
distances  decrease  significantly  when  the  flame  speeds  are  reduced  so 
much?"  Edse  emd  Lawrence  Jr.^  pointed  out  that:  after  a combustible 

gas  mixture  has  been  ignited  in  the  confines  of  a detonation  tube,  the 
greater  the  resistance  of  the  unburned  gas  to  the  motion  of  the  expanding 
burned  gas,  the  higher  will  be  the  pressure  of  the  burned  gas.  Thus  mod- 
erately strong  compression  waves  are  generated  in  the  \mburned  gas  which 
merge  rapidly  euid  thus  form  a detonation  wave  at  a very  short  distance 
from  the  ignitor. 
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Since  the  resistance  of  the  unbvtrned  gas  is  a function  of  its  density 
emd  friction  at  the  tube  wall,  the  greater  the  density  of  the  unburned 
gas  the  more  it  will  keep  the  burning  gas  from  expanding  so  that  pressxrre 
waves  are  formed  more  readily  and  thus  lead  to  a rapid  transition  from 
deflagration  to  detonation.  Since  the  lowering  of  the  initial  gas  tem- 
perature at  a constant  pressure  (l  atm)  increases  the  density  of  the  gas 
and  hence  the  resistance  to  the  flow  of  the  expanding  b\irned  gas,  rather 
strong  compression  vraves  are  generated  right  at  the  onset  of  ignition. 

This  fact  is  evidenced  by  the  pressure  traces  shown  in  Fig.  5«  At  123  K 
and  at  15  cm  from  the  ignitor  the  pressure  rise  is  already  of  the  order 
of  8 atm,  at  173  K it  amounts  to  5 atmospheres,  and  at  223  K it  is  only 
4 atm. 

As  shown  in  Fig.  7j  at  room  temperatvire  the  pressTore  is  essentially 
constant  even  up  to  150  cm  from  the  ignitor.  Thus,  these  strong  compres- 
sion waves  generated  at  low  temperature  increase  the  pressure  and  temper- 
ature of  the  unburned  gas  much  more  rapidly  than  at  room  temperature. 

In  addition,  the  reflection  of  the  precompression  shock  from  the  tube 
walls  and  the  downstream  end  of  the  tube  increase  the  temperatvire  and 
pressure  of  the  unburned  gas  even  further  by  increasing  the  strength  of 
the  precompression  shock.  An  increase  in  temperature  and  pressure  of 
the  compressed  but  unburned  gas  increases  the  bxirning  velocity  of  the 
ccanbustion  wave  very  rapidly.  Furthermore,  when  the  reflection  of  rather 
strong  compression  waves  passes  through  the  combustion  wave  a large  in- 
crease of  the  flame  surface  occurs,®  so  that  rather  turbulent  flames  are 
produced  which  also  lead  to  a significant  increase  of  the  burning  rate.® 
Also,  as  pointed  out  by  Brinkley  and  Lewis,''  the  effect  of  reflected  waves 
on  the  acceleration  of  the  flame  is  much  more  effective  than  the  incre- 
mental acceleration. 

According  to  Karlovitz,®  a very  strong  turbulent  wrinkling  of  the 
combustion  surface  can  be  produced  by  a very  strong  shock  wave  so  that 
a detonation  wave  is  formed  more  quickly.  Thus  as  long  as  the  burning 
speed  of  the  initial  gas  mixture  is  large  enough  to  maintain  a combustion 
wave,  the  much  higher  rise  in  pressure  and  temperature  of  the  compressed 
but  unburned  gas  due  to  the  lowering  of  its  temperature  along  with  the 
generation  of  turbulence  (as  explained  above),  is  always  going  to  increase 
the  burning  rates  to  a much  greater  extent  than  at  room  temperature. 
Consequently  a significant  reduction  of  the  induction  distance  is  accom- 
plished. 

To  explain  why  at  123  K a C.J.  detonation  wave  once  formed  is  unable 
to  propagate  further  down  the  tube  with  a constant  speed  and  pressvire, 
the  stability  and  structure  of  the  so-called  self- sustained  detonation 
wave  will  be  examined.  Many  papers  have  been  published  after  Zeldovitch,® 
von  Nevunan,^°  and  Doring^^  proposed  their  model  regarding  the  stability 
and  structure  of  self-sustaining  detonation  waves.  According  to  the  ZND 
model,  a gaseous  detonation  wave  consists  of  a shock  wave  which  traverses 
the  unreacted  gas  mixture  and  is  followed  by  a reaction  zone.  The  chem- 
ical changes  are  initiated  in  the  high-temperature  gas  behind  the  shock 
wave.  It  is  assvimed  that  at  the  tail  of  the  reaction  zone,  the  gas 
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reaches  thermodynamic  and  chemical  equilibrium.  The  energy  released  by  ; 

the  combustion  of  the  fuel  supports  the  propagation  of  the  shock  wave. 

Thus,  the  detonation  wave  (that  is,  the  combination  of  the  shock  and  com- 
bustion waves)  possesses  a self-sustaining  or  stable  character  and  prop- 
agates with  a constant  speed  and  pressure,  commonly  referred  to  as  the 
C.J.  velocity  and  C.J.  pressure.  The  detonation  wave  speed  and  pres- 
s\ire  are  independent  of  the  chemical  kinetics  of  the  combustion  reaction 
and  they  are  governed  only  by  the  fluid  mechanical  and  thermodynamical 
laws.  The  ZND  model  is  based  on  the  one-dimensional  structure  of  the 
detonation  wave. 

However,  as  observed  by  Fay^^  the  struct\ire  of  a self- sustained 
detonation  wave  is  strongly  influenced  by  the  chemical  kinetics  of  the 
combustion  process.  It  is  quite  certain^^  that  the  wave  front  of  all 
self-sustaining  detonation  waves  is  three-dimensional.  The  one-dimen- 
sional flow  behind  a plane  detonation  front  is  usually  disturbed  by  the 
occurrence  of  three-dimensional  oscillations  in  the  gas.  Martin  and 
White, using  interferometry  to  investigate  near  limit  detonations, 
foimd  convolutions  in  the  reacting  gases  which  they  attributed  to  spin. 

Dbring^^  has  suggested  that  at  the  limit  of  detonability  rarefaction 
waves,  generated  by  cooling  at  the  wall  of  the  tube,  may  penetrate  far 
enough  into  the  reaction  zone  to  quench  the  detonation  wave.  Sokolik^^ 
taking  into  consideration  Boring's  idea  that  a rarefaction  wave  is  formed 
at  the  tube  wall  dtiring  cooling  of  the  hot  gases  in  the  reaction  zone, 
postulated  the  following  criterion  for  the  decay  of  the  detonation  wave: 

d/2C  < ty  , 

where 

d = tube  diameter, 

c = speed  of  sound  (i.e.,  the  speed  with  which  the 
rarefaction  wave  will  travel),  and 


ri. 


1 

; t 


tj.  = reaction  time. 

Zeldovitch®  also  considers  the  heat  losses  from  the  hot  boundary 
layer  to  the  cold  wall  of  the  detonation  tube.  This  heat  loss  lowers 
the  temperature  in  the  reaction  zone,  which  causes  a decrease  of  the 
chemical  reaction  rates  and  thus  may  lead  to  an  incomplete  liberation  of 
heat  from  the  combustion  process.  Consequently  the  strength  of  the  shock 
wave  will  be  reduced,  which,  in  turn,  will  reduce  the  reaction  rate  still 
mere.  As  the  reaction  rate  is  decreased  the  thickness  of  the  reaction 
zone  is  increased.  The.  effect  of  thermal  losses  as  explained  by  Boring, 
Sokolik,^®  and  Zeldovitch®  applies  only  to  one-dimensional  detonation 
waves.  The  three-dimensional  analysis  of  the  wave  includes  heat  losses 
automatically. 
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Then  after  the  detonation  wave  has  decayed  to  a deflagration  wave, 
suitable  conditions  of  temperature  and  pressure  farther  down  the  tube 
may  cause  the  flame  to  accelerate  again  and  form  another  shock  wave  which 
results  in  the  formation  of  a second  stable  detonation  wave.  The  high 
temperature  and  pressure  may  be  caused  by  reflection  of  the  initial  shock 
wave  from  the  closed  end  of  the  tube.  This  explanation  is  in  agreement 
with  the  data  given  in  Table  I;  in  one  experiment  a c.J.  detonation  was 
observed  at  577  cm  frcra  the  ignitor  after  the  decay  of  the  first  wave 
which  was  formed  at  230  cm  from  the  ignitor.  It  appears  that  this  behavior 
of  the  detonation  wave  occurs  only  at  123  K because  of  the  higher  thermal 
losses  at  this  low  temperature. 


II.  CALCULATION  OF  THERMODYNAMIC  AI®  GAS-DYNAMIC  PROPERTIES 
OF  A ONE-DIMENSIONAL  GAS  FLOV/  TO  V7HICH  ONE  OR  MORE  DIF- 
FERENT GASES  ARE  ADDED 


A.  INTRODUCTION 

The  addition  of  gases  to  a primary  flow  through  a duct  may  be  ac- 
complished either  with  an  increase  in  duct  area  as  shown  in  Fig.  8 or  in 
a constant  area  duct  as  depicted  in  Fig.  9* 

In  the  following  analysis  it  will  be  assumed  that  the  flow  is  adia- 
batic relative  to  the  environment  and  that  there  are  no  viscous  inter- 
actions between  the  flow  and  the  walls  of  the  duct.  The  calculation  of 
the  length  of  the  channel  required  to  produce  a completely  homogeneous 
gas  mixture  will  not  be  included  in  this  analysis  because  this  problem 
involves  the  complex  nonequilibrium  and  rate  phenomena  of  diffusion,  heat 
transfer,  momentum  exchange,  and  chemical  changes.  Tiie  inclusion  of 
these  relationships  leads  to  very  complex,  and,  in  all  but  the  simplest, 
highly  idealized  cases  only  to  approximate,  equations  which  cannot  be 
solved  rigorously  at  the  present  time. 


B.  DERIVATION  OF  EQUATIONS 

For  the  derivation  of  the  equations  which  must  be  used  to  calculate 
the  thermodynamic  and  gasdynamic  properties  of  the  uniform  final  flow 
(subscript  m when  no  chemical  changes  occur  and  subscript  c when  chemical 
changes  have  occurred  and  produced  a mixture  in  chemical  equilibrium) 
from  the  initial  conditions  (subscript  i)  of  the  unmixed  gases,  the  fol- 
lowing notations  and  definitions  will  be  used: 

Ao  = area  of  duct  through  which  the  primary  flow  passes. 

Subscript  0 (zero)  is  used  to  denote  all  properties  of  the  primary 
flow 


FOR  NORMAL  INJECTION  ALL  a,=7T/2  AND  THUS  A, 


Addition  of  gases  to  a one-dimensional  flow  without  change  of  duct  area  (A, 


Aj  = area  of  injector  through  which  the  species  is  (i  = 1,  or  2, 

3,  n)  added  to  the  primary  flow  (i  - O)  at 

ai  = angle  between  direction  of  primary  flow  and  that  of  axis  of 

injector  carrying  species  i 

i = subscript  to  denote  chemical  species  of  gases,  primary  flow 

(i  = 0)  and  injected  gases,  i = 1,  2,  3 n or  i = 1,  2, 

3>  z when  chemical  changes  are  involved 

n = nxomber  of  injectors 

z = n\amber  of  species  in  mixture  in  chemical  equilibria 

Tn  1 = molecular  mass  (kg/kmol)  of  species  i (i  = 2,  3>  n 

or  z) 

Tn  „ = molecular  mass  of  gas  mixture  when  no  chemical  changes  are 

involved 

"Tn  g = molecular  mass  of  gas  mixture  when  chemical  changes  have  occur- 
red and  complete  chemical  equilibriian  has  been  established 

Tj  = static  temperature  of  species  i (i  = 0,  1,  2,  n)  just  prior 

to  mixing  (K) 

T„  = static  temperature  of  uniform  final  flow  when  no  chemical  changes 
have  occurred 

Tg  = static  temperatvire  of  \miform  final  flow  when  chemical  changes 
have  occurred  and  chemical  equilibrium  has  been  established 

Similarly  we  have: 

p = static  pressure  with  subscripts  i,  m,  or  c (N/m^) 

u = linear  gas  speed  with  subscripts  i,  m,  or  c (m/s) 


Tijn  *=  mole  fractions 


of  unreacted  gas  mixtvire 


Xij„  = mass  fractions) 

qi,c  = mole  fractions  of  reacted  gas  mixture 
Vi  = number  of  moles  of  species  i 
H = absolute  molar  enthalpy  (j/mol) 
h = absolute  enthalpy  of  a unit  mass  [J/kg  = (m/s)®] 


Eo,i  = molar  zero  point  energy  of  species  i 
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1 = formation  enthalpy  of  species  i at  temperature  T (j/kmol)  from 
* elements  at  T = OK. 

JR  = 8314.33  (j/kmol*K)  = universal  gas  constant  (j/kg*K) 

Ri  = (JRM)1 


= (9l/Z»m)  > specific  gas  constants  (j/kg»K) 

Rc  = («/^c)J 

The  three  fundamental,  unknown  properties  (namely,  temperature, 
speed,  and  pressure)  of  the  flow  of  the  mixture  can  be  calculated  from 
equations  derived  from  the  conservation  laws  of  mass,  momentvun,  and 
energy,  and  the  equation  of  state  in  the  forms  in  which  they  apply  to 
this  kind  of  flow. 

The  first  relationship  between  the  known  initial  and  the  unknown 
final  conditions  of  a flow  is  obtained  from  the  law  of  conservaticxi  of 
mass  (continuity  equation). 

For  the  mass  flow  rate  through  a channel  we  can  write 

m = puA  (1) 

When  p and  u represent  average  values  of  the  density  and  speed  over 
the  cross-sectional  area  A. 

Upon  elimination  of  the  density  by  means  of  the  equation  of  state 

P . . (2) 

we  obtain 

. p.7?7.u.A  / » 

® “ JR  * T ' 3 ) 

According  to  the  law  of  conservation  of  mass  we  have 


m = E “i  » 

i=0 


which  together  with  Eq.  (3)  leads  to 

n 

rmi 


p _ = Im . 

" A„  977  „u„  777  A„  u„977„A„  L.  T* 

i=0 
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or 


■aa 


n 

iricUcA^ 


= 


Uc 


z 

i=0 


771 


i^iAi 


(6) 


n 

where  A„  = Aq  + Ajcosai  or  A„  = Aq  for  the  case  of  normal  in- 
i=l 

jection,  or  as  shown  in  Fig.  9* 

The  second  relationship  between  known  and  unknown  variables  is 
based  on  the  conservation  of  mcanentum.  As  long  as  all  viscous  interac- 
tions between  the  fluid  and  the  walls  of  the  duct  and  injectors  through 
which  they  pass  can  be  disregarded,  the  total  mcanentvun  of  the  fluid,  or 
more  precisely  its  mcmentum  flux  at  any  cross  sectional  area  in  the  duct, 
is  made  up  of  two  contributions;  i.e.,  (l)  the  momentum  flux,  mu,  re- 
sulting from  the  motion  of  the  fluid  as  a whole,  and  (2)  the  momentum 
flux  caused  by  the  random  motion  of  the  individual  molecules  of  the  gas 
which  is  equal  to  the  product  of  the  static  pressure  and  the  channel 
cross-sectional  area.  Tius,  for  such  flows,  conservation  of  momentum 
leads  to 


n 

“o"o  + AqPo  + £ (miUi  * + At  pi)  = 

i=l  ' 

+A„p„  = BiUc  +A„pc  . (7) 

With  Aijx  = Aicostti,  Ui^x  = Ujcosoi  and  Eq.  (3)  the  momentum  equation 
can  be  written  in  the  following  form; 


The  areas  A©,  Aj,  and  A„  can  be  removed  frcsn  this  expression  by  using 
the  continuity  equation  in  the  form 

Substitution  leads  to 
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where  Xq  = (mo/jo)  and  the  Xi  = (mi/m)  are  the  so-called  mass  fractions 
of  the  individual  constituents  of  the  initial  flow.  Equation  (lO), 
which  satisfies  both  the  law  of  conservation  of  momentum  and  that  of 
mass,  contains  only  two  unknowns  (namely,  T„,  and  u„,)  when  no  chemical 
changes  occur,  since  in  this  caseW„  can  be  calculated  readily;  whereas 
the  final  gas  ccmposition  must  be  determined  first  to  calculate 

= E^i,c77?i>  where  the  qi  c are  the  mole  fractions  of  the  individual 
constituents  of  the  final  mijrtxu-e  in  chemical  equilibri\im.  Because  of 
the  greater  complexity  of  these  calculations  this  case  will  be  treated 
later . 

1.  Mixing  Without  Chemical  Changes 

Turning  now  t6  the  case  of  chemically  inert  mixtures,  we  find  that 
the  molec\ilar  mass  of  the  mixture,  '77?^,  can  be  calculated  readily  by 
means  of  the  known  mass  fraction,  xi,  of  the  individual  flows.  Fre- 
quently it  is  advantageous  to  use  the  mole  fractions,  for  these  cal- 
culations . 

By  definition  we  have 

(11) 


where  the  represent  the  number  of  moles  of  species  i passing  through 
channel  Aj  (i  = 0,  1,  2...n)  per  unit  time.  These  mole  numbers  are 
equal  to  the  mass  flow  rate,  mi,  divided  by  the  corresponding  molecular 
mass,  Wi,  of  the  constituent  according  to  the  equation 


Substitution  of  Eq.  (l2)  into  Eq.  (ll)  leads  to 


The  molecular  mass  of  a mixture,  by  definition,  is 

= I (14) 

or 

w. . f:  .uw, . (15) 

i=0 


Substitution  of  Eq.  (13)  into  Eq.  (15)  leads  to 


For  numerical  calc\ilations  Eq.  (lO)  is  written  as  an  equation  giving  the 
temperature  T„  as  a function  of  the  gas  speed  u„;  i.e.. 


t(M-C)  _ an_(  Um^i^  - Ur.  ) ^ (17) 

Rm 

Where 

'4^max  = *0^10  ^ 1 + ^ f 1 ^ ^ ' 

Alternately,  Eq.  (10)  can  also  be  written  as  an  equation  giving  the  gas 
speed  as  a function  of  the  temperature;  i.e.. 

The  superscript  M-C  is  used  to  indicate  that  numerical  values  of  the 
temperature  calculated  by  means  of  Eq.  (17)  for  certain  speeds,  u„,  or 
numerical  values  of  the  speed  calcvilated  by  means  of  Eq.  (18)  for  cer- 
tain temperatures,  T„,  satisfy  only  the  momentum  and  continuity  equa- 
tions. To  calculate  the  temperature  T„  Tor  the  gas  speed  u„) 
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which  satisfies  also  the  energy  equation  it  is  necessary  to  use  the 
energy  equation  to  remove  the  gas  speed,  u„,  frcan  Eq.  (17).  In  other 
words  we  must  obtain  an  equation  which  contains  T„  as  the  only  unknown 
or  we  must  remove  the  temperature  T„  frcan  Eq.  (l8)  to  obtain  an  equation 
which  contains  u„  as  the  only  variable . 


Because  of  the  complicated  relationship  between  the  thermodynamic 
functions  (such  as  enthalpies,  internal  energies,  and  specific  heats) 
and  the  temperature,  a rigorous  solution  can  be  established  only  for 
calorically  perfect  gases  [temperature-independent  specific  heats,  see 
Eq.  (37)].  For  all  other  gases,  even  when  no  chemical  changes  occur, 
the  final  conditions  (u„,  T„,  and  p„)  can  be  calculated  only  by  an  iter- 
ation procedure.  The  accuracy  of  the  resiilts  can  be  as  good  as  desired. 
Of  course,  it  cannot  exceed  that  of  the  thermodynamic  functions  used  in 
these  calculations. 


As  long  as  the  flow  is  adiabatic  relative  to  its  environment  the 
energy  equation  for  the  mixing  process  can  be  written  in  the  following 
simple  form: 


#) 


(19) 


The  injection  angles  Oi  do  not  appear  in  the  energy  equation  because 
they  do  not  have  any  influence  on  the  scalar  property  of  energy.  The 
terms  representing  the  heat  transfer  between  the  various  gas  flows,  which 
is  necessary  to  achieve  a xaiiform  temperature,  do  not  have  to  be  included 
in  the  energy  equation  because  only  the  final  state  of  the  gas  is  to  be 
calculated.  This  internal  transfer  of  energy  is  included  implicitly  in 
Eq.  (19).  It  affects  only  the  magnitudes  of  the  two  terms  representing 
the  energy  of  the  flow;  i.d.,  h,  the  enthalpy  which  consists  of  the  ran- 
dom motion  energies  of  the  molecules  plus  the  flow  work  term  pv,  and 
(u^/2),  the  kinetic  energy  due  to  the  motion  of  the  gas  as  a whole.  In 
the  same  manner  the  viscous  interactions  between  the  various  flows  are 
taken  into  consideration  in  the  momentum  equation  in  >rtiich  viscous  in- 
teractions determine  the  relative  magnitudes  of  the  momenta  due  to  the 
random  motion  (p)  and  ordered  motion  (pu^)  to  the  total  mcxnentum  which 
remains  constant  as  long  as  wall  friction  is  disregarded. 

In  Eq.  (19)  the  absolute  enthalpies,  hi,  of  species  i (usually  a 
single  component  gas;  if  not,  such  as  air,  it  may  be  treated  as  if  it 
were  a pure  gas)  are  related  to  the  absolute  molar  enthalpies  of  the 
species  by  the  following  expression: 


(20) 


>rtiere  the  superscript  Ti 
the  absolute  enthalpy  is 


denotes  the  temperature  at  which  the  value  of 
to  be  taken. 
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On  the  other  hand  the  absolute  enthalpy,  h^">,  of  a unit  mass  of  a 
mixture  is  obtained  by  summing  the  properly  weighted  molar  enthalpies, 
of  the  species;  e.g., 


h^ni  _ 53  n 


where  T„  is  the  temperature  of  the  gas  mixture,  which  means  that  all 
constituents  have  this  temperature.  The  weighting  factors  in  Eq.  (21), 
the  rii  m,  are  the  so-called  specific  molalities  of  the  constituents. 
They  are  defined  as  follows; 


ni,„  = 


moles  of  species  i 1 
.unit  mass  of  mixture J 


Using  Eq.  (12)  we  can  express  the  specific  molalities  in  terms  of  the 
mass  fractions; 


As  long  as  only  steady  flows  are  considered  it  is  not  necessary  to 
distinguish  between  mass  (mj  or  m)  and  mass  flow  rate  (mi  or  m)  in  many 
of  these  relationships.  Therefore,  to  simplify  the  notation,  the  dot 
over  the  properties  denoting  flow  rates  will  be  omitted  when  permissible. 

Sufficiently  accurate  numerical  values  for  the  absolute  molar  en- 
thalpies of  the  species,  Hi',  are  not  available  because  they  can  be  ob- 
tained only  from  determinations  of  the  molecular  masses  extending  to  at 
least  l6  significant  digits,  which  is  many  orders  of  magnitude  beyond 
our  present  capabilities.  The  absolute  enthalpies  include  the  rest  mass 
of  the  species,  expressed  in  units  of  energy,  Eo,i  = I’C®  (c  = speed 
of  light)  which  make  the  numerical  values  so  enormously  large.  Since  we 
are  never  interested  in  absolute  enthalpies,  because  only  differences 
appear  in  all  equations  involving  enthalpy  (or  internal  energy) , the  ab- 
solute values  can  always  be  replaced  by  properly  selected  relative  values 
of  enthalpy  (or  internal  energy) . When  no  chemical  changes  (which  means 
no  mass  changes  according  to  Einstein's  mass-energy  equivalence  relation- 
ship, Ae  = -An»c^)  occur  in  the  process  under  consideration,  the  absolute 
molar  enthalpies,  h|,  can  be  replaced  by  the  sensible  enthalpies, 

(H-Eq)T*  Ihus  using  the  relationships  Eqs.  (20)  and  (21)  and  introducing 
the  sensible  enthalpies,  the  energy  equation  [Eq.  (19)]  can  be  written 
as  follows; 

io  **  * ^.1 ' " " f]  ■ 

For  ccmputational  purposes  it  is  advantageous  to  replace  the  sensi- 
ble enthalpies  by  the  reduced  sensible  enthalpies,  ( » t>ecause 
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they  are  more  readily  available  in  tables  than  the  sensible  enthalpies. 
Therefore,  Eq.  (24)  is  modified  by  miiltiplying  the  first  term  in  each 

square  bracket  by  the  identities  or  respectively.  Furthermore, 

Jvl  Jvj.  ^ 

the  whole  equation  is  divided  by  the  mass  of  the  mixture,  m,  and  the  spe- 
cific molalities  are  eliminated  by  means  of  Eq.  (23)  to  express  the  com- 
position in  terms  of  the  mass  fractions.  We  arrive  at 


Solving  for  the  gas  speed  we  obtain 


The  superscript  (e)  is  used  to  indicate  that  this  gas  speed,  when  calcu- 
lated for  an  arbitrary  value  of  the  temperature  T„,,  satisfies  only  the 
energy  equation.  Usually  its  n\americal  value  is  obtained  with  an  arbi- 
trarily assumed  temperature,  T„.  The  last  term  under  the  square  root 
sign  in  Eq.  (25)  can  be  simplified  slightly  by  replacing  the  mass  frac- 
tion, Xi,  by  the  mole  fractions,  qi.  According  to  Eq.  (13)  we  have 


and  according  to  Eq.  (l4)  we  have 


(26) 


Upon  division  of  Eq.  (26)  by  Eq.  (27)  we  obtain 


(27) 


An  approxjjnate  value  of  the  gas  speed,  u„,  can  be  calculated  when  the 
expression  for  the  temperatvire  of  the  gas  mixture,  according  to 

Eq.  (17),  is  substituted  into  Eq.  (29).  The  resulting  quadratic  equation 
in  u„  may  be  considered  to  contain  u„  as  the  only  unknown  because  the 

term  V qj  ( . ^ ) varies  only  very  slightly  with  temperature.  For 

i^  \ /i 

calorically  perfect  gases  the  reduced  sensible  enthalpies  are  constant 
so  that  for  such  gases  the  speed  of  the  gas  mixtures  can  be  calculated 
directly  from  the  given  conditions.  For  gases  with  temperature-dependent 
reduced  sensible  enthalpies  the  speed  of  the  mixture  can  be  calculated 
by  an  iteration  procedure  starting  with  a reasonable  estimate  of  T„  so 

n /H-E 

that  the  sum  V qj  ( ~ ^ ) can  be  determined.  This  procedure  will  be 
i=0  ' ^ 

developed  after  a general  insight  into  the  relationships  between  T„  and 
u„  has  been  gained. 


For  an  efficient  application  of  the  iteration  procedure  and  for  a 
critical  appraisal  of  the  solution  obtained  by  the  method  described  above 
it  is  helpful  to  examine  the  possible  range  of  values  of  T„  and  u„  as 
given  by  the  combined  momentiam  and  continuity  equations  alone  [Eqs.  (17) 
or  (18)]  and  by  the  ene-’gy  equation  alone  [Eq.  (29)]. 

According  to  the  momentum  equation  [Eq.  (17)]  the  range  of  possible 
speeds  of  the  gas  mixture  is  given  by  the  following  condition; 


0 S u„ 


^ (M-C) 

2 nnijinax 


V u 1 


cosai 


For  gas  speeds  outside  of  this  range  the  temperature  of  the  gas  mixture 
would  be  negative,  an  obviously  impossible  result.  Therefore,  the  term 
in  the  square  bracket  of  Eq.  (30)  represents  the  maximvun  possible  speed 
of  the  gas  mixture  according  to  the  momentum  equation.  It  is  obvious 
that  T,^””^)  becomes  zero  when  u„  = 0 and  when  u„  = v^,max*  The  maximum 
speed,  and  the  specific  gas  constant,  R„,  are  constants  for  a 

given  set  of  initial  conditions  (xj,  uj,  Ti,  Wj,  and  tti).  An  inspection 
of  Eq.  (17)  [or  (I8) ] reveals  at  once  that  the  relationship  between  u„ 
and  T„  satisfying  the  raomentvun  and  continuity  equations  simultaneously 
is  that  of  a parabola  and  that  the  temperatvire  of  the  gas  mixture  attains 
the  largest  possible  value  for  a given  set  of  initial  conditions  when 
u„  = 1/2  Thus  we  have 


^(M-C)  (o.^Um|.max)  . 

^m,max  n 

Furthermore,  we  note  [see  Eq.  (I8)]  that  for  each  possible  temperature 
of  the  gas  mixture  there  are  two  different  speeds  which  the  mixture  may 
attain.  This  fact  requires  special  attention;  It  Is  necessary  to  know 


29 


whether  both  solutions  are  physically  possible  and  what  processes  they 
represent. 


Turning  nov?  to  the  energy  equation  [Eq.  (29)]  we  recognize  at 
once  that  the  range  of  possible  speeds  of  the  gas  mixture  is  given  by 
the  following  condition: 


0 S u„  s 


The  maximum  gas  speed,  u„,  according  to  the  energy  equation  is  attained 
when  the  temperature  of  the  mixture  becomes  zero.  Thus  we  write 

II  X,  [uf  . 2H.T,  . 4%.  ■ (33) 

The  gas  speed  becomes  zero  when  the  temperature  reaches  its  highest  pos- 
sible value,  which,  according  to  Eq.  (29),  is 


t(e)  _ 
■‘^m,max  - 


(um!^max)^ 

1=0 


/h-Eq 

V iRT  /i 


The  relationship  between  temperature,  T„,,  and  speed,  u„,  according  to 
the  energy  equation  can  now  be  written  in  the  simple  form 


,(E)  _ 


- 2H.T.  i ,,  . 

i=0  ^ 


According  to  this  equation  the  speed  of  the  mixture  decreases  steadily 
as  the  temperature  is  increased  from  T„  = 0 to  its  maximum  value  as 
given  by  Eq.  (34). 

For  the  case  of  a simple. tangential  (aj  = 0)  injection  of  hydrogen 
into  a flow  of  oxygen  the  u^®^  versus  T„  [Eq.  (35)1  relationship,  to- 
gether with  the  versus  T„  [Eq.  (l8)]  relationship,  are  depicted 

in  Fig.  10  for  two  different  injection  speeds  (uu  = 200  m/s  and 

2 

Ujj^  = 700  m/s).  Since  the  actual  flow  of  the  gas  mixture  must  satisfy 

all  three  fundamental  equations  (namely,  continuity,  momentum  [Eqs.  (l?) 
or  (l8)],  and  energy  [Eq.  (35)1)  only  those  points  common  to  both  cxirves 
represent  possible  conditions  of  the  flow  of  the  gas  mixture. 

Figure  10  reveals  that  there  are  4 possibilities;  i.e., 

(l)  the  two  curves  do  not  intersect  at  all  as  in  the  case  for  which 
= 700  m/s. 
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O =200  m/s  0.1836 

A = "^00  m/s  =0.6426  ( no  solution ) 

^Og  ■ I^Og  =0.8735 


PHYSICALLY 

IMPOSSIBLE 

SOLUTION 


T„  =366K 
= I354m/s 
=2.385 


\ \ 


Z' 


T„=7I6K 
=422  m/s 
M„  =0.5305 


Fig.  10  - Relationships  between  final  temperature,  T„,  and  gas  speed 
u„,  of  mixture  (subsonic  flow) 


versus 


r 


i 

1. 


(2) 

(3) 

(4) 


the  versus  T„  curve  is  tangent  to  the 

T„  curve  (uj^  = 4^.5  m/s), 

the  two  curves  intersect  at  two  points  (uj^  = 200  m/s),  and 
the  two  curves  intersect  only  at  one  point  '^^inax^' 


When  the  two  ciirves  have  no  point  in  common,  there  is  no  solution  because 
the  mixing  of  the  flows  leads  to  choking.  In  other  words  the  mixture 
cannot  pass  through  the  channel  without  proper  adjustment  of  the  initial 
conditions . 


When  the  two  curves  have  only  the  ixjint  in  common  where  they  are 
tangent  to  each  other  (possibility  No.  2),  the  Mach  number  of  the  mixture 
is  one  (see  Example  No.  1).  The  point  of.tangency  occurs  always  at  a 
temperature  which  is  slightly  less  than  Tjj^^j^ax  ^ speed  which  is 

somewhat  greater  than  1/2  Um,max-  For  the  case  that  the  final  gas  tem- 
perature is  equal  to  T^jinax  the  Mach  number  of  the  flow  is  less  than  one, 
as  the  following  calculation  shows; 


j^im,max  _ 


1 „(M-C) 

2 ^m,max 


1 „(M-C) 

2 *^,max 


I (M-C)  I ^ (o.^  J7 


(36) 


where  the  value  for  given  by  Eq.  (31)  has  been  used.  The  point 

of  tangency  may  be  considered  as  the  critical  point  since  it  is  the  lim- 
iting point  of  unchoked  flow.  Choking  would  occur  at  once  when  the  speed 
of  one  of  the  two  gases  is  increased  slightly.  The  condition  of  choking 
can  also  be  produced  or  relieved  readily  by  variations  of  the  tempera- 
tures, injection  speeds,  mole  fractions  .of  the  constituents,  and  injec- 
tion angle.  It  is  obvious  that  the  choking  condition  is  more  quickly 
produced  with  normal  injection  because  of  the  constant  area  condition 
(see  Example  No.  3)- 


When  the  two  curves  have  two  points  in  common,  there  seem  to  be  two 
different  final  flows  for  a given  mixing  process.  For  initially  subsonic 
flows  the  solution  with  the  higher  temperature  T„,  and  lower  speed  u„ 

(see  Fig.  10)  represents  a subsonic  flow  (Mn,  < 1),  whereas  the  solution 
with  the  lower  temperatxire  T„  but  higher  speed  u„  represents  a supersonic 
flow  (M„,  > l) . To  determine  whether  both  solutions  represent  physically 
possible  processes  it  is  necessary  to  calculate  the  entropy  changes  re- 
sulting from  these  processes.  These  calculations  reveal  that  steady 
subsonic  flows  cannot  be  made  supersonic  by  mixing  because  this  transi- 
tion involves  a spontaneous  decrease  in  entropy.  Such  a process  could 
be  accomplished  only  by  an  expansion  shock  wave  which  is  physically  im- 
possible. On  the  other  hand  an  initially  supersonic  flow  may  remain 
supersonic  or  it  may  become  subsonic  by  experiencing  a normal  shock 
trajisition  (see  Example  No.  2). 


32 


I 

d 


I 


When  the  initial  speeds  are  lov/  or,  more  generally  speaking,  when 
the  ordered  motion  kinetic  energy  and  momentum  are  negligible  in  compar- 
ison with  those  of  the  random  motion,  the  maximum  possible  speed  from 
the  momentvim  equation  [Eq.  (30)]  becomes  larger  than  that  derived  from 
the  energy  equation  [Eq.  (33)1  and  consequently  the  two  curves  have  only 
one  point  in  common  in  the  positive  coordinate  quadrant.  In  this  case 
only  the  subsonic  solution  exists;  the  second  solution  (supersonic; 
expansion  shock!)  involves  a negative  temperature. 

In  many  respects  the  u^®^  - T„  and  - T„  curves  may  be  con- 

sidered as  the  counterparts  of  the  Rayleigh  and  Hugoniot  curves  (in  the 
p-v  plane)  for  constant  area  duct  flows  with  heat  exchange  or  the  Fanno 
curves  for  constant  area  duct  flows  subjected  to  viscous  interaction  at 
the  duct  walls.  However,  it  is  impossible  to  state  how  the  various 
points  of  the  u„  - T„  curves  can  be  obtained  by  a certain  mixing  process 
because  of  the  large  number  of  initial  parameters  (uj,  Ti,  Xj,  'TTTj,  and 
«i)  whereas  Rayleigh,  Hugoniot,  and  Fanno  curves  result  from  the  varia- 
tion of  a single  parameter  (heat  exchange  or  friction).  Therefore,  any 
change  of  any  of  the  parameters  of  the  initial  gas  (xi,  Ui,  Ti,  W^) 
changes  both  the  ui^)  - T„  and  the  - T„  curve  (see  Fig.  lO)  whereas 

a change  of  q does  not  affect  the  Rayleigh  line. 

According  to  the  previous  discussions  it  is  obvious  that  the  tem- 
peratiire  T„  and  speed  u„  of  the  gas  mixture  can  be  calculated  by  means 
of  several  different  iteration  procedures.  The  basic  problem  consists 
in  finding  a solution  for  two  variables  which  appear  in  two  equations 
which  must  be  satisfied  simultaneously  but  cannot  be  combined  readily 
because  of  the  complicated  mathematical  relationships  between  these  two 
variables . 


The  solution  requiring  the  lowest  number  of  iterations  is  based  on 
the  fact  that  a rather  reasonable  value  of  the  gas  speed,  u„  can  be  cal- 

culated  when  it  is 


assumed  that  the  term  ^ tjj  y ) ~ changes 


only  slightly  with  temperature  and  that  the  temperat\ire  of  the  mixture, 
T„,  cannot  exceed  the  lower  of  the  maxima  given  by  Eqs.  (30)  and  (33). 
Substitution  of  t^^"^)  from  Eq.  (17)  into  Eq.  (35)  leads  to  a quadratic 
equation  in  u„  with  the  solution 


(M-C)  E 
Um  = Uni, max  ‘ 


1 ± 


/ (E) 

_ |Um.niax 

IuCm-c) 

Y^,max 


where 


(37) 


^ " i?o  ^ 


(38) 
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For  subsonic  flows  the  ix)sitive  sign  in  front  of  the  radicand  leads  to 
the  physically  impossible  expansion  shock  wave.  Furthermore,  we  note 
that  Eq.  (37)  has  real  solutions  only  when 


•Ibis  condition  may  serve  as  a simple  critericaa  to  determine  whether 
choking  occurs  or  not.  For  a quick  answer  to  this  question  Fig.  11  was 
prepared.  It  depicts  the  dimensionless 
u /u^^’ 

speed  ^ function  of  I 1 various  values  of  E.  For  a 


speed  as  a function  of  I C)  ) various  values  of  E.  For  a 

^m,max  \'^,max/ 

calorically  perfect  gas  E is  independent  of  temperatvure  and  no  iteration 
is  required.  When  E depends  on  temperature,  the  accuracy  of  the  gas 
speed  calculated  by  means  of  Eq.  (37)  depends  on  the  accuracy  of  the 
estimated  temperatTire  and  the  rate  at  which  E depends  on  temperature. 
With  the  fairly  close  approximation  to  the  correct  gas  speed  of  the  mix- 
t\ire,  as  given  by  Eq.  (37),  the  temperature,  as  given  by  the  momentum 
equation  [Eq.  (17)],  is  calculated.  This  temperature  is  now  used  as 


the  starting  point  of  the  iteration  procedure  by  calculating  the  differ- 
ence u,\^]  - = AU-.  by  means  of  Eq.  (35)  and  Eq.  (l8) . For  subsonic 

flows  and  for  all  u„  < i/2  'this  difference  is  positive  when  the 

temperat\u'e  is  less  than  the  correct  value  and  vice  versa.  For  subsonic 
flows  and  u„  > 1/2  Uj^^“2x  “this  difference  is  positive  when  the  tempera- 
ture, used  for  calculating  this  difference,  is  greater  than  the  correct 
value  and  vice  versa  (see  Fig.  10  and  Table  6).  The  calculations  are 
repeated  with  properly  adjusted  temperatures  t1“) , until 

L(e)  . u(M-C)|  g 

where  e is  a small  number  which  should  not  be  smaller  than  warranted  by 
the  acc\iracy  of  the  thermodynamic  functions  used  in  these  calculations. 

After  the  correct  values  of  the  final  gas  speed,  u„,  and  temperature, 
T„,  have  been  determined,  the  static  pressure  of  the  gas  mixtvire,  p„,  can 
be  calculated  by  means  of  Eq.  (5)  from  which  the  channel  cross-sectional 
area,  A„,  is  removed  by  the  following  relationships: 


Am  = Ao  + Ajcosai  > 
i=l 


* - p,7n,ui  ’ 


Table  6.  Iteration  for  calculation  of  T„  and  u„ 


(a)  = 492  m/s 

2 

with  = 650  K we  obtain  from  Eq.  (37)  = 736.8  m/s; 

with  K we  obtain  from  Eq.  (37)  ul^^  = 7^5  ta/s. 

Then  from  Eq.  (I7):  t1^^  = ^ ( 1303.. 8 ^u^  = 643  K 

■^m 

and  64l  K,  respectively. 


ijt(n) 

*^in 

A 

A u 

643 

737.47 

738.30 

0.83 

642 

741.18 

741.52 

0.34 

64l 

744.75 

744.73 

-0.02 

640 

748.19 

747.92 

-0.27 

639 

751.50 

751.10 

-o.4o 

638 

754.71 

754.26 

-0.45 

637 

757.82 

757.41 

-0.4l 

636 

760.85 

760.55 

-0.30 

635 

763.79 

763.67 

-0.12 

634 

766.65 

766.78 

0.13 

t-  subsonic 
solution 


-supersonic 

solution 


Subsonic  Solution:  T„  = 641.05  K,  u„  = 744.56  m/s,  M„  = 0.9840 

Supersonic  Solution:  T„  = 634.48  K,  u„  = 765.24  m/s,  M„  = I.OI6 
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Table  6.  (Continued) 


Ujj  = 100  m/s 

with  700  K we  obtain  from  Eq.  (37)  ui^) 

Then  from  Eq.  (l7)  T^^^  = 735-7  K 


= 292.3  m/s. 


JE) 

Au(") 

(K) 

(m/s) 

(m/s) 

(m/s) 

735 

291.990 

284.266 

-7.724 

734 

291.507 

292.628 

1.121 

734.125 

291.567 

291.596 

0.029 

Solution:  T„  = 734.13  K,  u„  = 291.58  m/s,  M„  = O.36II 


so  that 


• 9^0  . • OTi 


'^n.U„ 


Upon  division  of  numerator  and  denominator  by  m and  cancellation 
of  91  we  arrive  at 


Po'^o'lo  ^ 

According  to  Eq.  (28)  we  have 

iil-  = JU-  . 

777 1 777„ 

Substitution  of  this  relationship  into  Eq.  (4l)  leads  to 
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^ ^qTo  + V iiZl  cosai 
L^oUo  PlUi 


The  ratio  of  the  duct  area  through  which  the  final  flow  of  the  mixture 
(ih,A„)  passes  to  that  through  which  only  the  primary  flow  (mo>Ao) 
passes  can  be  calculated  as  follows: 


from 

Am  = Ao  + 23  AiCOSOi  , 

i=l 

with 

A _ .MTq. 
PoUoWo 

and 

A,  - , 

^ PiUiWi 

we  obtain 

cosa, 

Ao  mo  To  Pi  Ui  ^i 


With  Ai  = m*Xi  and  Xi  =•  rn 
we  arrive  at 


The  increase  in  entropy  caused  by  the  mixing  of  different  gas  flows 
follows  frcm  the  expression  for  the  entropy  of  gas  mixtures.  Before 
mixing  we  have 


s = 23  xiSi  » 

i=0 


where 


After  mixing  we  have 
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_ - Y.  Tie’ll  - ^Pn 

1-0  \ 

/I  i=0 

and  thus 


o o ” /oP=i\'^n>  /oP=i\^i  n n 

“ 1?0  J ■ S ■ i?0  ^ 

2.  Mixing  of  Gases  at  Rest 

The  previous  equations  can  be  used  also  to  calculate  the  temperature 
resvilting  from  the  mixing  of  gases  at  rest  (uq  = 0 and  all  ui  = 0).  When 
the  gases  are  not  flowing,  only  the  energy  equation  is  available  for  cal- 
culating the  final  temperature  T„.  When  the  pressure  is  held  constant, 
Eq*  (33)  can  be  used.  For  Uq  = Ui  = 0 we  obtain 

- ^2  t (!^)"‘  . (48) 

Substitution  of  this  expression  into  Eq.  (34)  leads  to 

E x,R  T ( 

,pstatic,dp=0  _ ^ ^ . (49) 

R r „ fH-Eor'" 

R„  Z-  1U  I gjrp 


^ nn^ 

as  given  by  Eq.  (28),  we  obtain 


T,static,dp=0  _ i=0  V £RT 

X — rk  / V 


At  constant  volume  the  internal  energies,  Ej,  instead  of  the 
enthalpies  must  be  used'.  With  E = H - IRT  we  obtain 


i 

i 
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,pstatic,dv=0  _ J-=o 

' I,  •'  (^): 


n 

. 1 

V t«r  /i 

V s 

RT  /i 

" X 

i=0 

Lv 

- 

I.  '■  («■)’ 


3.  Mixing  of  Gas  Flows  Involving  Chemical  Changes 

When  chemical  reactions  occur  as  the  gases  are  mixed,  the  calcula- 
tions become  much  more  tedious.  First  of  all  the  reduced  sensible  en- 
thalpies in  the  energy  equation  must  be  replaced  by  the  reduced  formation 


enthalpie 


Only  the  case  that  the  final  gas  mixture  is  in  a 


state  of  complete  thermodynamic  and  chemical  equilibrium  will  be  dis- 
cussed here.  Since  the  composition  of  the  reacting  gas  mixture  depends 
not  only  on  temperature  but  also  on  pressure,  the  calculations  must  be 
started  with  estimated  values  of  the  final  temperature,  T^,,  and  the 
final  pressure,  p^.  The  subscript  c is  used  to  stress  the  fact  that 
chemical  reactions  have  been  included  and  to  distinguish  this  state  frran 
the  non-reacting  mixtoire  for  which  the  subscript  m is  being  used. 

With  the  estimated  values  of  temperature,  and  pressvire,  p^^^, 

the  composition  (^i,c)  of  the  final  gas  is  calculated  according  to  the 
standard  procedure.  Then,  according  to  the  energy  equation  [see  Eq. 

(35) ]>  the  gas  speed  Uq,  is  calculated  for  this  temperature  and  composi- 


SRgTcZqi  c 


«m, limit  = ^ Xi  ['^1  + 2RiTi  J . (53) 

(e) 

For  exothermic  chemical  changes  the  speed  vte,^limit  is  no  longer  the 
largest  possible  value.  At  low  temperatures  the  sum  of  the  formation 
enthalpies  at  T^  in  Eq.  (52)  may  become  negative.  Therefore,  the  soeed 


'im,^limit  is  reached  at  a temperature  where  rii^c  is  equal  to 

fE^ 

zero.  For  lower  temperatures  u^  ’ exceeds  the  so-called  limiting  value, 
M,  limit*  Whether  these  conditions  represent  physically  possible  solu- 
tions cannot  be  said  in  general.  Since  they  can  occur  only  in  subsonic 
flows  (these  temperatures  are  lower  than  the  constant  pressure  adiabatic 
flame  temperatures)  it  may  be  expected  that  these  low  temperatures  lead 
to  the  impossible  expanslcn  shock. 


Since  the  calcxilations  can  be  started  only  with  an  arbitrarily 
assumed  pressure,  the  correct  pressure  must  be  established  by  iteration. 
According  to  the  continuity  equation  the  pressure  must  satisfy  the  condi- 
tion [see  Eq.  (3)] 

Pc  - A„  * uc-^c  ’ 

» 

where  ^ is  given  by  the  thermodynamic  and  gasdynamic  values  of  the  ini- 
tial  flovfs  before  mixing.  With 


we  obtain 


= Ao  + Ajcosaj 

i=l 


, mi  KTi 


„^o  + V ^ — 

Po  * ^0^0  ^2.  Pi  * ^1^1 


COSOi 


Upon  division  of  numerator  and  denominator  by  m and  with  Eq.  (28)  we 
arrive  at 


R«  ^ ^ cosai 

PoUo  Pi"l  _ 


where  S is  the  molecular  mass  of  the  unreacted  mixture. 

i=0  * 

Substitution  of  Eq.  (56)  into  Eq.  (5^)  leads  to 


^ ^ |_PoUo  PxUi 

Ulus  the  pressure  corresponding  to  the  assumed  temperature,  is 
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(58) 


P 


c,calc. 


For  the  selected  temperatxire  of  the  chemically  reacting  gas  mix- 

ture an  increase  in  pressure  reduces  the  amount  of  dissociation  as  long 
as  at  least  one  of  the  mole  number  changes,  Av(^',  of  the  cr  chemical 
changes  occurring  in  the  mixture  is  not  zero.  On  the  other  hand  a de- 
crease in  pressure  increases  the  amount  of  dissociation  in  this  case. 

z 


Because  of  this  effect  the  molecular  mass,  77^ 


= E 

i=l 


^i,c 


(z  = number 


of  species  in  equilibrium  mixture),  and  the  speed,  u^®),  will  increase 
as  the  pressure  is  increased  and  vice  versa.  Consequently,  the  pressure 
calculated  by  means  of  Eq.  (58)  will  decrease  when  a higher  value  is 

Therefore,  a new  estimate,  p^’^)  can 


used  for  the  estimated  pressure, 
be  obtained  as  a 


laabcu  cxux  e j ci  ucw  co  uxiucibc  j jj' 

simple  average  of  the  calculated  and  estimated ’values; 


4”Lt 


= ( 


,(n-l)  . „(n- 


Pc,caic  + Pc 


n-1)) 

,est/ 


0.5 


(59) 


When  the  cranposition  of  the  reacting  gas  mixture  is  not  affected 
very  much  by  pressure  variations,  the  iteration  formula  [Eq.  (59)]  may 
be  replaced  by  the  expression 


(60) 


where  x is  an  empirical  factor  between  0 and  0.5.  Tiie  iterations  are 
continued  \intil  the  difference  between  calculated  and  estimated  pressure 
satisfies  the  condition 


•>c”Lic  - 


s 6 


(61) 


where  c is  a small  number  in  accordance  with  the  desired  accuracy.  After 
the  correct  pressure,  p^,  has  been  determined.it  is  necessary  to  check 
the  validity  of  the  temperature. estimate,  To  perform  this  test 

both  gas  speeds,  u^®]  and  u(^“^),  are  calculated  with  this  temperature 
estimate,  Tl^',  and  the  corresponding  pressure,  p^,  by  means  of  Eq.  (52) 
and  Eq.  (l8),  which  now  is  written  as 


u 


(M-C)  _ 1 (M-C)  ^ 

0 - 2 ^ 


^c"^c  * 


(62) 


where,  as  in  the  non-reacting  case. 


I • ; 


•i  I 

i I 

I > { 


'4”m2x  = xqUo  (l  + xiu,  (l  + cosai  . 

However,  the  maxiraiun  temperatvire,  Tc^ma^j  inust  now  be  based  on  the  molec- 
\ilar  mass  of  the  mixture  of  the  products  formed  by  the  chemical  changes 
at  this  temperature.  Therefore,  instead  of  Eq.  (3I)  we  have  for  the 
reacting  mixture 

t(m-c)  ^ (M_c)  . (63) 

v-,mcxA  -^Cjmax 

Because  of  this  relationship  the  temperature,  can  only  be  esti- 

mated. Actually  a rigorous  value  is  not  needed.’  However  this  tempera- 
ture limit  may  be  of  some  help  in  establishing  the  first  estimate, 
which  must  be  lower  than  this  limit  although  it  is  not  possible  to  state 

how  much  lower  it  should  be.  If  = u^^^  - / 0,  the  calcu- 


lations have  to  be  repeated  with  T, 


(n+1) 


, when  ■'  > 0 for 


subsonic  flow  and  u^,  < 1/2  whereas  when  Au^”^^  > 0 

for  subsonic  flow  and  Ug  > l/2  u^^^ax*  When  < 0,  the  opposite 

changes  in  temperature  must  be  ma^e.  For  supersonic  flows,  positive 

values  of  Au^”^  indicate  that  the  temperatiire  estimate,  was  too 

low  and  vice  versa.  The  iterations  are  continued  until  |au^^'|  S € 
where  e is  again  a properly  chosen  small  value. 


TO  illustrate  these  techniques  let  us  consider  the  flow  of  oxygen 
(primary  flow)  to  which  hydrogen  is  added  at  various  conditions  (Example 
No.  1). 

>'pproximate  values  of  the  final  speed  of  the  gas  mixture,  u^^^,  as 
calculated  from  Eq.  (37)  with  estimated  temperatures  obtained  from 
Table  7>  are  compiled  in  Table  9 together  with  the  rigorous  calculations 

Am 

of  T„,  u„,  p„,  M„,  1 — , and  the  increase  in  entropy,  — . This  table 

reveals  that  approximate  values  of  the  gas  speeds  deviate  in  most  cases 
by  less  than  1^  from  the  exact  values.  The  estimated  temperatvires,  of 
course,  are  less  reliable.  Table  10  a and  b have  been  included  to  show 
the  rate  of  convergence  of  the  iterations  leading  to  the  correct  values 
of  u„  for  = 492  (m/s)  and  u^^  = 100  (m/s),  respectively.  For 

Ujjg  = 492  (m/s)  the  iteration  producing  the  supersonic  solution  (impos- 
sible expansion  shock)  has  been  included  in  order  to  show  that  this  in- 
jection speed  produces  two  solutions  which  are  very  close  to  the  critical 
conditions  (M„  = 1)  at  which  the  flow  becomes  choked.  For  Uj^  = 493  (m/s) 
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EXAMPLE  1 


A subsonic  flow  of  hydrogen  is  injected  tangentially  into  a subsonic 
flow  of  oxygen. 

The  initial  conditions  of  the  oxygen  flow  are: 

= 1400  K 
Vs 

M02  = 0.8735 

Uq^  = 600  m/s 

p-  = 1 atm 

O2 


= 0.36  (xq^  = 0.9)  = 32  kg/kmol 


The  initial  conditions  of  the  various  hydrogen  flows  are: 


- 

200 

K (for 

all  speeds  Ujj^ 

"Ha  = 

100 

m/s 

Mjj^  = 0.0918 

= 

200 

m/s 

= 0.1836 

= 

0 

0 

m/s 

= 0.3672 

= 

450 

m/s 

= 0.4131 

= 

490 

m/s 

= 0.4498 

= 

492 

m/s 

= 0.4517 

1 atm 


tiH„  = 0.64  (xjj„  = 0.1)  ^ Ho  = 2 kg/kmol 

/ (e) 

The  values  of  u^^niax>  '^jmax*  '^m,max>  ( ?M-C^  ) seven 

\uffl,max/ 

different  hydrogen  injection  speeds,  are  listed  in  Table  7*  For 

the  temperattires,  which  the  final  flows  are  estimated  to  attain,  the 

values  of  53  = &i,i  * functions  of  E appearing  in 

E<1*  (37)  are  compiled  in  Table  8. 
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with  = 640  K we  obtain  from  Eq.  (37)  =*  7^5  n>/s 

Then  from  Eq.  (17):  (j-303i8  « 643 


741.18 


744.75 


738.30 


741.52 


744.73 


subsonic 

solution 


748.19  747.92 


751-50 


751.10 


754.71  754.26 


757.82 


760.85 


760.55 


-0.41 


-0.30 


763.79  763.67 


766.65  766.78 


4-supersonic 

solutl(Hi 


Subsonic  Solution:  T.  - 641.05  K,  u«  - 744.56  >/s,  ■ 

Supersonic  Solution:  T.  ■ 634.46  K,  u.  ■ 765.24  n/s, 


Table  10.  (continued) 

UHg  = 100  m/s  ; 

with  = 700  K we  obtain  from  Eq.  (37),  = 292.3  m/s 

Then  from  Eq.  (17):  t1^^  = 735-7  K 


m 

(K) 

^M-C) 

(m/s) 

(e) 

(m/s) 

(n) 

Au 

(m/s) 

735 

291.990 

284.266 

-7.724 

734 

291.507 

292.628 

1.121 

734.125 

291.567 

291.596 

0.029 

Solution:  T„  = 73^.13  K,  u„  = 291.58  m/s,  M„  = O.36II 


there  is  no  solution.  Hence,  choking  begins  at  injection  speeds  between 
492  and  493  m/s.  Because  of  the  large  number  of  initial  parameters  of 
the  flows  it  is  not  possible  to  prove  analytically  that  at  the  point  of 
choking  the  two  cxirves  become  tangent  to  each  other.  However,  a graph- 
ical presentation  of  the  conditions  near  choking  in  the  u„  - T„  plane 
shows  that  the  condition  of  tangency  exists  at  M„,  = 1.  At  the  point 
where  the  temperature  tJ^M-C)  reaches  its  maximum  (ujj^  = 450  m/s)  we  have 

Mn  = ^ (see  Table  9)«  The  effect  of  the  hydrogen  injection  speed, 

U|T  , on  the  final  state  of  the  mixture  is  shown  in  Fig.  12. 
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EXAMPLE  2 


A supersonic  flow  of  hydrogen  is  injected  tangentially  into  a 
supersonic  flow  of  oxygen.  The  initial  conditions  are; 


2000  m/s 

^2 

= 1400  m/i 

II 

as 

200  K 

^0. 

= 1000  K 

1 atm 

POa 

= 1 atm 

= 

1.8287 

' 

M02 

= 2.3969 

II 

CM 

0.64 

^02, 

= 0.36 

: conditions  we  have 

- 

%i,max  - 

1 / 

= 23  XlUl  (1  + 

i=0  ' 

! 1668.60  m/s  , 

II 

= J E 'if 

V i=0  L 

+ 2RjTi 

(Si^aY 

V ST  /, 

rr 

1 

2121.17  m^ 

(M-c)  _ - 

m mov  “ ^ “ 


,(E) 


^m,max.  - 


- 1071.59  K , 

,2 

= 962  K,  and 


( 


V^m,max/ 


,^)Biax/ 


According  to  Fig.  11  solutions  exist  since  L > 3.5.  Assuming 

T„  = 600  K (z  = 3.5108)  for  the  supersonic  solution > we  obtain  from 

Eq.  (37)  (+  sign). 


= 1419  m/s  . 


With  this  value  the  first  temperature  estimate  according  to  Eq.  (17)  is 

= 545  K . 


F ^ 

\i 


i 1 


i ^ 
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;,-  ♦ 


K'l 
[;•  ^ 


t, 
pi 


^ 1 


Iteraticjn  leads  to  u„  = l4l7.72  m/s  and  T„  = 547.63  K,  M„  = 2.0224  so 
that  according  to  Eq.  (57)  P„  = 1.202814  atm.  Assianing  T„  = 900  K 
(Z  = 3.5803)  for  the  subsonic  solution  we  obtain  from  Eq.  (37)  (-  sign), 

ui^)  = 511  m/s  . 

With  this  value  the  first  temperature  estimate  according  to  Eq.  (17) 
becomes 

= 911  K . 

Iteration  leads  to  u„,  = 5IO.9I  m/s  and  T„  = 910.58  K so  that  K.  = 0.5701 
and  p„  = 5.549776  atm. 

The  ratio  of  the  pressures  of  the  two  solutions  is 

p(M»<l)  , ^ 

According  to  the  normal  shock  relationships  for  a gas  with  variable 
specific  heats  we  have 

tei..  V ftl./  

#1,.  te’-  - «■>)  ■ (E*  ■ » '))■  • (tt.,  ■ 

.mere  . 

With  Ta  = = 910.58  K and  Ti  = = 547.63  K , 

we  have  = 3.583  and  = 3.499  , 
so  that  (Ss\  = 4.6l4  , 

VVn.s. 

which  is  an  exact  agreement  with  the  ratio  of  the  pressures  of  the  two 
solutions; 

, , , 

fei)  ' 

^e  curves  of  versus  T„  and  u^®^  versus  T„  are  shown  in  Fig.  I3. 

Although  these  curves  represent  relationships  between  the  final 
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Fig.  13  . Relationships  between  final  temperature,  T.,  and  gas  speed, 
Ua,  of  mixture  (initial  flows  are  supersonic) 
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temperatvire  and  gas  speed,  only  the  points  ccxnmon  to  both  curves  are 
related  to  the  initial  conditions;  all  other  points  have  no  relationship 
to  them.  Also,  because  of  the  reduction  in  number  of  variables  from 
Ti>  Pi,  ui  to  T„,  p„,  u„  it  is  not  possible  to  present  grajdiically  the 
path  fran  the  initial  to  the  final  conditions. 


EXAMPLE  3 


A subsonic  flow  of  hydrogen  is  injected  normally  into  a subsonic 
flow  of  oxygen.  Ihe  initial  conditions  are 


Ujjg  = 200  m/s 

UQg  = 300  m/s 

%2  = 200  K 

= l400  K 

Pj^  = 1 atm 

p*  = 1 atm 

Mfig  = 0.182872 

1^2  = 0.436736 

TI02  " 

The  linear  speed  of  the  primary  flow  (oxygen)  in  this  example  had  to  be 
much  lower  than  that  in  the  first  example  (tangential  injection)  in  order 
to  avoid  choking  even  when  the  hydrogen  is  added  at  a very  low  speed. 

With  the  initial  conditions  as  given  above  we  have 

= ^2  • «0a  • " ^361.255812  m/s 


and 


^,max 


Xl 


uf  + 2R 


= 1795.22  m/s 


and  thus 


= 1.739B  . 


According  to  Fig.  11  this  speed  ratio  is  well  below  the  critical  value 
leading  to  choking  when  Z % 3*5. 


For  the  meocimum  temperatxire  we  obtain 


= 713.185023  K 


and 

= 702  K (Z  = 3-53376  @ 700  K)  . 

In  accordance  with  these  values  we  use  an  estimated  temperature  of 
600  K for  which  Z = 3-5108  and  thus  frcan  Eq.  (37) 

= 485-95  m/s  , 
and  according  to  Eq.  (l?) 

= 654.8  K . 


Iteration  leads  to 


and 


T« 


653.16  K 
483-19  m/s 


M„  = 0.632804 


r 


P»  = T.  ■ = 0.804623  atm  . 

U . 

" P02% 

An  increase  of  the  oxygen  speed,  Up^,  from  3OO  to  400  m/s,  reduces 
^^mai  ■4°  1178.44  m/s  and  increases  u^^^ax  1812.68  m/s,  so  that 

= 2.366, 


which,  according  to  Fig.  11,  means  that  the  addition  of  hydrogen  to  this 
flow  causes  choking.  On  the  other  hand,  an  increase  of  the  hydrogen 
speed,  Ujjg,  from  200  to  800  m/s,  raises  the  final  Mach  number  by  2.56^ 

(from  0.6328  to  0.6490).  Ihe  final  conditions  for  this  case  are 

T,  » 662.5  K 1 

I M,  = 0.6490 
u,  « 499  m/s  I 


Pa  ■ 0.7903  stm. 


EXAMPLE  4 


A subsonic  flow  of  hydrogen  = 0.0918)  is  added  tangentially  to 
a subsonic  flow  of  oxygen  a 0.2912) . Calculate  the  final  state  of 
the  mixture  after  combustion  has  occurred  and  complete  chemical  equilib- 
rium has  been  established. 

Because  of  the  high  combustion  enthalpy  of  hydrogen,  choking  occurs 
at  the  oxygen  speeds,  Uq^,  used  in  the  previous  examples.  Using  the 
criterion 


hij  (”-  ■ kJ 


which  applies  to  constant  area  duct  flows  with  heat  addition,  it  was 
estimated  that  the  following  initial  parameters  would  not  lead  to  choking; 


Ujj  = 100  m/s 


Uq  = 200  m/s 


Th^  = 200  K 
= 0.64 


Tog  = 1400  K 
Tio^  = 0.36 


The  final  state  of  the  gas  mixt\ire  was  calculated  by  two  methods: 

(l)  First  the  state  of  the  mixtiure  was  calculated  without  consid- 
eration of  any  chemical  changes  according  to  the  procedure  explained  in 
the  previous  examples.  Starting  then  with  these  uniform  conditions 
(denoted  by  the  subscript  m)  the  final  state  (denoted  by  the  subscript 
c)  was  determined  by  finding  for  a constant  area  duct  flow  (A^  = Am)  the 
pressxire  and  temperature  ratios  which  satisfy  simultaneously  the  Rayleigh 
line  (momentiun  equation  from  which  u^.  has  been  removed  by  means  of  the 
continuity  equation), 


Bi-1 


1 - Ici 
v_ 


= (s.r.)‘ 


and  the  Hugoniot  curve  (energy  equation  frcm  which  both  speeds  have  been 
eliminated  by  means  of  the  mcxnentum  and  continuity  equations), 


t 
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where 


.2m 

V„  T„  p^. 


(67) 


Both  the  Rayleigh  and  Hugoniot  relationships  contain  only  Tg  and  pg 
as  the  \onknown  variables.  However,  an  explicit  solution  cannot  be  es- 
tablished because  of  the  complicated  relationships  between  these  variables 
and  the  formation  enthalpies 


_T_ 

T„ 


7T? 


(68) 


Therefore,  the  two  equations  (65  and  66)  must  be  solved  simultaneously 
by  an  iteration  technique  which  has  to  be  started  not  only  with  an  esti- 
mated vsLlue  of  the  final  temperature,  Tg,  but  also  with  an  estimate  of 
the  final  pressure,  p . 


With  these  estimates  first  the  composition  of  the  final  gas 

is  calculated  (see  Worksheet  No.  1).  Using  this  cca&position  we  must 
determine  whether  the  estimated  pressure,  together  with  the  estimated 
temperature,  satisfies  the  Hugoniot  equation  (66).  This  test  is  per- 
formed by  determining  a so-called  calculated  presstire  fran  Eq.  (66) . 

In  deriving  the  expression  for  the  calculated  pressure,  Pgjgaic* 
feet  of  pressvire  on  the  formation  enthalpy  of  the  final  gas  mixture  is 
ignored.  Thus  from  Eq.  (66)  we  obtain 


where 


and 


4%ic  - ■> 


a + 7e,est  . ^ 


a « - 0. 

R»T„ 


, PJc 


‘] 


(69) 


(70) 


The  superscript  (n)  denotes  the  nth  estimate.  If  p(>^)  , ^ p(n)  the 

calculation  of  the  composition  is  repeated  with  adjusted  values  of  the 
estimated  pressure,  but  the  same  temperature  t(^^ . Since  the  ef- 

fect of  pressure  on  the  formation  enthalpy  of  the  final  gas  mixtvire  is 
usually  not  very  large,  the  calcvilated  pressure  is,  in  most  cases,  fairly 
close  to  the  correct  presstire  so  that  successive  estimates  may  be  obtained 
by  the  formula 


p(n) 

'^c,est 


+ X 


(n) 

c,calc 


5(n)  . 

c,est 


) 


56 


WORKSHEET  NO.  1 

Procedure  for  Calctilatlog  the  Mole  Fractions,  Tji,  of  a Combustion  Gas  in  Equilibrium 
Containing  the  Species:  HsO,  Ha,  Oa,  OH,  0,  and  H 


.6  atm  P,  » 1.0042  atm 


.40). 

PJ'A  = 0.145185 

T^^>  - 3000  K 

K(HaO) 

• P^  « 17.063590 

Fq  • 0.6  atm 

K 

- 1.1967 

v§ 

- 0.5625 
''Ha 

. k(h)  . 

PJ'A  » 0.203203 

M.  - 0.2298 

Assume 


' 


[d]  + [I)*]  + 1 - E] 


- 0 - m 


• □ ' 5oa 

- 0.078400 

32.000 

• ® “ *V) 

> 0.040652 

-»  d)  16.000 

13  • ® • 4<XI 

- O.U6158 

□1  17.008 

• ® -nji 

- 0.070443 

0 1.008 

• 0 ■’'h* 

- 0.120175 

□]  2.016 

S ' © " ’hfeO 

= 0.574172 

til  18.016 

7r\.  . 2:n»,7W, 

- 1.000000 
■ 15.792409  kg/kmol 

^ ([2'*'  0)'O’5  + S'*"  ^ 

If  Q]  ^ 0,  Repeat  Calculations  with  - EJ  • 2)  -♦  0 

If  More  Than  Two  Iterations  are  Required  to  Make  I 0 1 < 0.001,  Use  Linear  Interpolation 
as  Follows  (For  n k 2): 


LiiLki  fil. 


where  0.5  S X < 1.  The  iterations  are  repeated  until 

l^c.calc  where  e is  again  a small  number  which  must  be 

chosen  properly.  A rather  accurate  value  of  the  pressure  can  be  obtained 
by  linear  interpolation  between  two  pressures  for  which  the  differences 

P^“calc  - P^^lst  = Pc^caic  " ^c^U  ^ 

case  we  have 


Pc  ■ P&sl  - 


(n-1) 


To  begin  the  calcvilations  with  temperature  and  pressure  estimates 
which  guarantee  a rapid  convergence  to  the  correct  values,  we  calculate 
first  the  adiabatic  flame  temperature  (see  Table  U)  of  the  mixture 
(initial  state  m)  at  constant  pressure  when  the  flow  is  subsonic 
(M„  < 1)  and  at  constant  volume  when  the  flow  is  supersonic  (M„  > l) • 

For  subsonic  flow  this  temperature  applies  to  the  case  that  s.r.  = 0 
[see  Eq.  (65)]  and  for  supersonic  flow  the  constant  volume  = ^)  adi- 

abatic  flame  temperature  makes  the  speed  ratio  (s.r.)  infinitely  large. 
For  subsonic  flows  the  adiabatic  flame  temperat\ire  at  constant  press\ire 
is  only  slightly  larger  than  T„,  and  p„  is  only  slightly  less  than  p„ 
because  heat  addition  to  a subsonic  flow  has  a mildly  expansive  effect, 
whereas  for  supersonic  flows  T^  is  much  higher  than  the  adiabatic  flame 
temperature  at  constant  volume,  and  p^  is  many  times  that  of  p„  because 
heat  addition  to  a supersonic  flow  has  a strongly  compressive  effect. 

After  a ccmpatible  pair  of  ^ and  ^ values  has  been  obtained  frcm 

Eq.  (69)  by  iteration,  it  is  used  to  calculate  the  dimensionless  speed 
of  the  initial  gas,  s.r.,  according  to  Eq.  (65)*  If  this  calculated 
speed  ratio  (or  its  square)  differs  frcm  the  actual  speed  ratio  of  the 


initial  gas 


all  previous  calculations  must  be  repeated  with  an 


improved  estimate,  T^  of  the  final  temperature,  together  with  a new 
value  for  p^.  For  subsonic  flows  the  speed  ratio,  s.r.,  is  increased  as 
Tg  ' is  decreased  and  vice  versa.  Therefore,  we  have 


est  > 


(s.r.)2<^,  (72) 

and  vice  versa.  Since  there  are  two  solutions,  it  is  necessary  to  cal- 
culate the  Mach  number.  Mg,  of  the  final  flow  to  make  sure  that  it  is 
less  than  one.  The  solution  with  Mg  > 1 has  no  physical  significance 
because  the  transition  from  subsonic  to  supersonic  flow  by  means  of  heat 
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Table  11.  Calculation  of  adiabatic  flame  temperatiire 


Tijj  = 0.64  = 0.36  p„  :=  1.0042  atm  T„  = 685.2  K 


= 12.81024  kg/kmol 


= 8.063347  • For  the  case 


that  u„  = 0 and  p^  = p„  we  have  Tg  = 3115«9  K as  shown  by 
the  table  below: 


= 3115.9 


addition  in  a constant  area  duct,  although  mathematically  possible,  can 
occur  only  when  an  expansion  shock  wave  is  involved  which  entails,  how- 
ever, a spontaneous  decrease  in  entropy. 

For  supersonic  flows  both  solutions  are  physically  possible.  One 
solution  represents  supersonic  flows  (M^  > l)  for  which  the  temperatures 
and  pressures  are  only  moderately  higher  than  those  of  the  adiabatic 
combustion  at  constant  volume.  For  this  solution  the  speed  ratio  is 
Increased  as  is  decreased,  and  vice  versa,  so  that  the  relationships 

[Eq.  (71)  and  (72)]  apply  here  too.  For  the  second  solution  leading  to 
Mc<l,  because  of  a normal  shock  wave,  the  final  temperature  and  pres- 
sure are  much  higher  than  those  of  the  adiabatic  flame  temperature  at 
constant  volume.  For  this  solution  the  speed  ratio,  s.r.,  is  increased 
as  the  temperature,  is  increased,  and  vice  versa.  Therefore,  it- 

erations for  this  solution  must  be  based  on  the  conditions 


t(“=1)  > * 

■‘c,est  ■^c,est  » 
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and  vice  versa.  The  iterations  are  continued  until 


(s.r.)2  - 


u: 


R «T 

IT 


where  again  e is  a small  number  chosen  in  accordance  with  the  accuracy 
desired  and  the  limitations  given  by  the  thermodynamic  data  (formation 
enthalpies  and  equilibrium  constants)  used  in  the  calculations.  A rath- 
er accurate  value  of  the  final  temperature  can  be  obtained  by  linear 
interpolation  between  two  temperatures  for  which  the  differences 
u® 

(s.r.)^  - s 6 are  small.  In  this  case  we  can  write 


Tc  = Tc^est  - 


(n)  (n-l) 


(2)  In  the  second  method,  the  details  of  which  are  given  in  Sec- 
tion 3>  the  calculation  of  the  state  of  the  unreacted  mixture  is  by- 
passed. Therefore,  it  is  more  difficult  to  begin  the  calculations  with 
a reasonable  estimate  of  the  final  temperature  and  pressure.  Consequently 
a larger  n\miber  of  calculations  is  required,  as  evidenced  by  the  data 
shown  in  Fig.  14.  The  results  of  both  techniques  are.  tabulated  in  Table 
12.  An  increase  of  the  speed  of  the  oxygen  flow  from  200  m/s  to  222  m/s 
produces  choking.  The  final  conditions  for  u = 221  m/s  are  given  in 
Table  13.  The  physically  impossible  solution  involving  the  subsonic- 
supersonic  transition  is  included  to  show  how  close  this  flow  is  to  that 
at  which  the  Mach  number  of  the  combustion  gas  is  one  (M^  = 1) . 


III.  TRANSITION  FROM  DEFLAGRATION  TO  DETONATION 


Most  previous  research^®"^®  on  three-dimensional  detonation  waves 
dealt  with  the  decay  of  waves  generated  by  ignitors  having  high  power 
densities.  Since  these  observations  are  not  necessarily  pertinent  to 
the  prediction  of  the  transition  from  deflagration  to  detonation,  exper- 
iments were  carried  out  in  the  Aeronautical  and  Astronautical  Research 
Laboratories  of  The  Ohio  State  University  to  study  the  propagation  of 
deflagrations  in  practically  unconfined  hydrogen-air  mixtures.  Hydrogen- 
air  mixtures,  contained  between  0.004  inch  polyethylene  sheets  draped 
over  a metal  frame,  were  ignited  in  the  center  by  an  electrically  heated 
copper  wire  having  a diameter  of  0.003  inch.  The  energy  transferred 
from  this  igniter  to  the  combustible  gas  mixture  was  sufficiently  low  to 
allow  the  combustion  process  to  start  as  a deflagration.  Furthermore, 
no  metal  vapor  or  hot  particles  were  propelled  through  the  combustible 
gas  mixtures,  which  in  all  experiments  contained  29.5^&  hydrogen.  The 
flames  were  photographed  on  Kodak  35ram  Tri-X  film  by  a FASTAX  camera. 
According  to  Figs.  15  and  l6  the  flame  speeds  are  given  by  the  following 
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Relationships  between  final  temperature,  Tg,  and  gas  speed,  Ug,  of  burned  gas  In  equilibrium; 
Initial  flow  Is  subsonic  (Mt,  = 0.09l8>  = 0.2912) 


Table  12.  Mixing  of  a subsonic  flow  of  oxygen  with  a subsonic  flow  of 
hydrogen  followed  by  combustion 


Initial  State  of 
Gases  before  Mixing 

Mixing  without 
Chemical  Change 

Combustion  Gas  in 
Complete  Equilibrivim 

Uq^  = 200  m/s 

= 179.56  m/s 

Uj.  = 957.05  m/s 
(953) 

Ujj^  = 100  m/s 

Tq^  = 1400  K 

T„  = 685.2  K 

T„  = 3044.5  K 
(3046) 

Tjjg  = 200  K 

P-  = 1 atm 
Oa 

p„  = 1.0042  atm 

p = 0.6881  atm 

® (0.690) 

p„  = 1 atm 
Ha 

Mog  = 0.2912 

KU  = 0.2298 

Mg  = 0.6801 

Mjjg  = 0.0918 

^Ha  " 1*^386 

7„  = 1.3721 

B 

7^  = 1.2972 

Oa 

^ Oa  = 32  kg/kmol 
^ Ha  = 2 kg/kmol 

„ = 12.8  kg/kmol 

c = 15.5718  kg/kmol 

1^,1  * 

TV,  a 0.64 

■Hajm 

"ib.c  ' 

The  values  in  parentheses  were  obtained  by  direct  calculation  (see 
Section  3) • 
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Table  13.  Mixing  of  subsonic  flow  of  oxygen  with  a subsonic  flow  of 
hydrogen  followed  by  combustion 


Initial  State  of 
Gases  before  Mixing 

Mixing  without 
Chemical  Change 

Combustion  Gas  in  Complete 
Chemical  Equilibrium 

U02  - 221  n/e 

u.  ■ 191.20  a/s 

„ , 1282  m/s 

Ug^  - 100  a/s 

(1390) 

Tq^  - UWX)  K 

T,  • 686.0  % 

- a 2995  K 

® (2976)  K 

Tgg  - 200  K 

POa  - 1 

Pa  ■ 1.0076  ata 

B a 0.535  ata 

- 1 ata 

(0.489)  ata 

0.3217 

M.  • 0.24U5 

« a 0.9205 
(1.0028) 

>%a  - 0.0918 
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Fig.  15  - Geometrical  representation  of  how  an  image  is  recorded  on  film 


Film  Plone 

Direction 
of  Film 
movement 


Fig.  16  - Trace  of  flame  on  rotating  film  strip 


equation: 


77?uJ  ■ dt 


tana  , 


where  b is  the  distance  between  the  camera  lens  and  the  ignitor,  F is 
the  focal  length  of  the  lens,  Tu  the  temperature  of  the  vmbxirned  gas, 

Tb  that  of  the  burned  gas,  and'^b  corresponding  molecular  mas- 
ses, ^ the  film  travel  per  unit  time,  and  a the  half-angle  of  the  cone 

Clu 

formed  by  the  flame  images  on  the  film  (see  Fig.  17). 

Measurements  were  made  in  bags  containing  1,  2,  and  4 m®  of  the 
hydrogen-air  mixture;  the  results  are  shown  in  Fig.  l8.  Whereas  the 
flame  propagation  rates  in  the  1 and  2 m®  bags  were  quite  reproducible, 
a fair  amount  of  scatter  prevailed  in  the  4 m^  volume.  The  most  signif- 
icant observation,  however,  is  the  fact  that  there  is  a noticeable  ac- 
celeration of  the  flame  propagation  (see  also  Fig.  17).  This  accelera- 
tion, although  quite  modest  in  the  1 m^  bag,  is  large  enough  in  the  4 m^ 
bag  that  it  may  be  assvuned  that  a greater  flame  speed  increase,  ixjssibly 
resulting  in  a detonation  wave,  would  occur  in  much  larger  volumes  of 
the  combustible  gas  mixture.  Unfortvinately,  experiments  with  larger 
bags  cannot  be  carried  out  at  this  laboratory  because  of  the  strong  com- 
pression waves  which  accompany  these  explosions  and  whose  strengths  in- 
crease rapidly  w-th  bag  size. 


Experiments  to  measure  the  pressure  rise  at  various  distances  from 
the  bags  are  in  progress. 


IV.  QUENCHING  OF  DETONATION  AND  DEFLAGRATION  WAVES 


Because  of  the  low  gas  speeds  and  small  pressure  changes  produced 
by  deflagration  waves,  the  principal  function  of  a flame  arrestor  con- 
sists of  terminating  the  progress  of  the  chemical  reactions  in  the  flame. 
On  the  other  hand  a detonation  wave  is  followed  by  the  high-speed  flow 
of  the  hot  combustion  gas  which  is  much  more  difficult  to  attenuate. 
Furthermore,  this  wave  involves  a large  rise  in  pressure  which  must  be 
dissipated  by  the  arrestor  in  order  to  prevent  ignition  of  the  downstream 
gas  by  the  shock  compression.  Therefore,  arrestors  designed  to  quench 
a detonation  wave  must  not  only  stop  the  progress  of  the  combustion  proc- 
ess but  they  must  also  produce  a drastic  decay  of  the  shock  wave  and 
prevent  the  hot  gases,  which  traverse  the  arrestor,  from  igniting  the 
downstream  combustible  gas  mixture. 

To  study  both  effects  independently,  two  different  experimental 
Investigations  have  been  initiated.  First  a shock  tube  has  been  designed 
cmd  constructed  for  measurements  of  the  effectiveness  of  arrestors  in 
bringing  about  a significant  decay  of  the  shock  wave  so  that  the  wave 
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emerging  fran  the  arrestor  can  no  longer  ignite  the  downstream  gas. 
Progress  on  this  investigation  has  been  slow  because  several  graduate 
students,  after  working  on  this  problem,  have  decided  to  discontinue 
their  studies  toward  an  advanced  degree. 

Another  apparatus  (see  schematic  views  shown  in  Figs.  19  and  20) 
has  been  designed  and  constructed  for  the  study  of  flame  quenching  by 
arrestors  made  of  different  materials,  having  various  porosities  and 
configurations.  The  test  section,  as  shown  in  Fig.  20,  accommodates  the 
arrestor  and  has  a rectangular  cross-section.  It  is  equipped  with  large 
windows  to  photographically  record  the  processes  occurring  at  the  ig- 
nitor surfaces.  Eight  photo  transistors  (Texas  Instrument  LS  400)  will 
be  used  to  determine  the  flame  propagation  rates.  Severed  thermocouples 
will  be  placed  in  the  upstream  and  downstream  sections  for  temperature 
measurements.  The  pressures  in  both  sections  will  be  monitored  by  quartz 
crystal  transducers.  Photographic  records  will  be  obtained  with  a 35nim 
Fastax  camera.  Hie  combustion  will  be  initiated  by  two  spark  plugs  fed 
by  a 7.5  hV  ignition  transformer.  It  is  planned  to  use  hydrogen-air, 
methane-air,  and  propane-air  mixtures  for  these  experiments.  Several 
arrestor  materials  have  been  given  to  us  by  the  Air  Force  Aero  Propulsion 
Laboratory.  These  samples  include  polyurethane  foams  of  different  pore 
sizes,  and  mesh  and  sponges  made  of  aluminum  and  stainless  steel. 

In  both  types  of  experiments  the  initial  gas  mixtures  are  at  rest. 
Future  experiments  are  planned  with  flowing  gas  mixtures. 


V.  PERFORMANCE  OF  SUPERSONIC  RAMJET 


Previous  calculations  of  the  thirust  specific  fuel  consumption  of 
supersonic  ramjets  (see  AFOSR-TR  76  O503)  revealed  that  the  best  per- 
formance can  be  attained  only  with  extremely  high  pressures  in  the  can- 
bustion  chamber  unless  operation  at  altitudes  where  p„  < lO""*  atm  ia 
considered. 

To  investigate  the  effects  of  (l)  less  efficient  diffusion  prv  ces- 
ses, (2)  the  use  of  hydrocarbon  fuel  (propane)  instead  of  hydrogen  suid 
(3)  the  mixing  of  fuel-air  streams  at  different  injection  speeds  and 
temperatures  on  the  thrust  specific  fuel  consvunption  of  the  engine  ad- 
ditional calculations  were  made,  the  results  of  which  are  compiled  in 
Table  l4.  According  to  these  data  highly  efficient  diffusers  are  not 
necessarily  desirable  when  high  flight  speeds  are  considered  because 
they  lead  to  very  high  pressures  in  the  combustion  chamber.  Although 
an  inefficient  diffusion  process  (e.g.,  a normal  shock  at  the  diffuser 
inlet)  can  be  used  to  avoid  the  high  static  pressure  in  the  combustion 
chamber,  the  losses  caused  by  the  rise  in  entropy  of  the  air  lead  to 
very  high  values  of  the  thrust  specific  fuel  consvunption.  An  addil'onal 
increase  in  the  thrust  specific  fuel  consumptlcxi  is  produced  when  tue 
temperature  and  the  injection  speed  of  the  fuel  differ  greatly  from  the 
respective  values  of  the  air  flow  through  the  combustion  chamber  inlet. 
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19  - General  arrangement  of  the  15  cm  I.D.  combustion  tube  to  be  used  for  flame  quenching 


GLASS  FRAME  *46.5  cm 

—TOTAL  GLASS  LENGTH  *42.5  cm 

5LASS  SECTION  FOR  PHOTOGRAPHY =40cm 
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